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Spin Fluctuation-Induced Superconductivity 
in k-BEDT-TTF Compounds 
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Spin fluctuation-induced superconductivity in quasi-two dimensional organic compounds, k- 
BEDT-TTF salts, is investigated within a fluctuation exchange (FLEX) approximation using a 
half-filled Hubbard model with a right-angled isosceles triangular lattice (transfer matrices — t, 
— t'), extending a previous work above T c . An energy gap of A2 or (x 2 — y 2 )-type develops with 
decreasing temperature below T c more rapidly than in the BCS model. The calculated dynam- 
ical susceptibilities enough below T c show sharp resonance peaks like those in certain cuprates 
superconductors. The calculated nuclear spin- lattice relaxation rate 1/Ti shows a T 3 behavior 
below T c in accordance with experiment. Estimated values of 1/Ti are roughly consistent with 
experimental results. A prediction is made for the doping concentration dependence of T c and 
the antiferromagnetic and superconductive instability points are calculated in the U /r vs. t'/t 
plane. 

KEYWORDS: organic superconductor, spin fluctuation, Hubbard model, triangular lattice, resonance peak, nuclear 
spin-lattice relaxation rate 
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Recent theoretical investigations indicate that the su- 
perconductivity in quasi-two dimensional organic com- 
pounds k-(BEDT-TTF) 2 X [where X = Cu{N(CN) 2 }X', 
X' = CI, Br] can hepftxpJained in terms of the spin fluctu- 
ation mechanismBa'&QEP The calculated values for the 
transition temperature T c as well as the d-wave charac- 
ter of the order parameter seem-to be-cansistent with the 
existing experimental results.BB§BII9'IId) 

In a previous paper (hereafter referred to as I) we have 
reported the results of calculations for T c and the normal 
state properties of these compounds by using a Hubbard 
model consisting of dimer orbitals.Ef The purpose of the 
present communication is to report on the results of our 
extended calculations mainly to include the supercon- 
ducting state. 

The model and the approach are the same as those dis- 
cussed in I. We use a half-filled Hubbard model with a 
right-angled isosceles triangular lattice consisting of the 
anti-bonding dimer orbitals with the inter-dimer trans- 
fer integrals — r and — t' as shown in Fig. 1. The esti- 
mated values for_th£ transfer integrals are r ~ 0.07eV 
and t'/t ~ 0.8.E3E3 We use the fluctuation exchange 
(FLEX) approximation where the dynamical suscepti- 
bilities are calculated within the renormalized random 
phase approximation (RRPA) and the normal and the 
anomalous self-energies are approximated by the sim- 
plest ones including a spin fluctuation propagator with- 
out vertex corrections. This approximation has been 
successful in|j^aja^pxe»oaj^Bi2Kstigations on bigh-T c 
cuprates, etcH-BBBBBBEl) and recent studies on 
the vertex corrections seem to indicate that, its effects 
will not seriously modify the FLEX results.E3 
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Fig. 1. (a) The model unit cell and the transfer integrals, (b) 
Unperturbed Fermi surface for r'/r = 0.8. Dashed lines show 
the antiferromagnetic zone boundary. 



Before reporting on the results below T c we show in 
Fig. 2 the calculated phase boundaries, in the U/t vs. 
t'/t plane, between the paramagnetic metallic (PM) and 
the antiferromagnetic (AF) phases and between PM and 
the superconductive (SC) phases. The former is obtained 
from the points where the value of I /xq extrapolated to 
T = vanishes. The latter is obtained from the points 
of vanishing T c calculated by extrapolating the T c /t vs. 
U/t plots, shossin in the inset of Fig. 2, using the Pade 
approximantsEiJ' This result, being consistent with the 
conjectures in I, shows that the AF and SC instabilities 
compete with each other and the latter wins for the val- 
ues of t'/t larger than ~ 0.3. In other words, the d-wave 
superconductivity appears near the antiferromagnetic in- 
stability only when the electronic structure is favorable; 
an antiferromagnetic instability is not necessarily accom- 
panied by a d-wave superconductivity. With increasing 
U /t from a superconducting instability point we expect 
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Fig. 2. Instability points of the paramagnetic state against super- 
conductivity and antifcrromagnetism. Inset shows the results of 
an improved calculation for T c /r vs. U/r. 
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Fig. 3. Electron number (doping concentration) dependence of 
Tc/r for t'/t = 0.8. 



a phase transition to an antiferromagnetic state. The 
phase boundary between the superconducting and anti- 
ferromagnetic phases must be obtained by comparing the 
free energies of the two phases, since the phase transition 
is naturally of the first order. This task is left for future 
studies. 

Figure 3 shows the calculated values of T c for t'/t = 
0.8 for various values of n, the number of electrons per 
site, and U /r. These results may be regarded as a pre- 
diction for the doping concentration-dependence of T c in 
k-(BEDT-TTF) 2 X. 

Now we discuss the properties in the superconducting 
state. We first show in Figs. 4(a) and 4(b) the anisotropy 
and the temperature dependences of the gap function 
Z\fc = Rc{Z\(fc, respectively, where the gap func- 
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Fig. 4. Calculated gap function, (a) Temperature dependence of 
Afc, with k indicated by an open circle on the Fermi line shown 
in the inset, (b) Anisotropy or the wave vector-dependence of 
Afc, with k on the Fermi line shown in the inset of (a). 



tion is defined by 
A(k,uj) = 



Z(k,Lj) : 



cvZ(k, lo) = w - - S {1) {k,uj + iri) 



irj) 



(1) 



XA 1 ) and S^> being the normal and the anomalous self- 
energies, respectively. The symmetry of the gap function 
is clearly of A 2 or (x 2 — y 2 )-type and its amplitude devel- 
ops more rapidly than in the standard BCS model below 
T c - Figure 5 shows the renormalized density of states 
at various temperatures, indicating how it is influenced 
by the d-wave gap formation with lowering temperature. 
As is expected from Fig. 4 the energy gap in the density 
of states develops rapidly with decreasing temperature, 
approaching close enough to a limiting result at around 
T c /2. 

The calculated dynamical susceptibilities well below T c 
show strong resonance peaks around (tt, it) and (tt, —it), 
just as in the calculations for the high-T c cuprates. Sim- 
ilarly to those in cupcates. we interpret this resonance 
peak as a spin excitonJISB a bound pair of an electron 
and a hole excited across the energy gap. Figure 6 shows 
the imaginary part of the dynamical spin susceptibility 
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Fig. 5. Density of states at various temperatures. 
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Fig. 6. Calculated imaginary part of the dynamical susceptibility 
for q = (tx, 7r) at various temperatures. 



at q = (ir, 7r) at various temperatures. We see how the 
resonance peak develops from a broad spectrum in the 
normal state. Figure 7 shows the dispersion and broad- 
ening of the spin excitons. The resonance peak position 
ojq and the width j q are given by 



l-URex s (q,u q ) = 0, 
{dRex s (q,w)/duj} 



(2) 



with 



X s (<Z, i^m) = -j-r ^2 I G ( fc + 1' ' luJ n + iw m ) G (fc, iw„) 



N, 



-F (k + q, \uj n + ioj m ) F (fc, iu n )] 



(3) 



where uj m and oj n are the Bose and Fermi Matsubara fre- 
quencies, respectively, and G(fe,ia; n ) and F(k,ko n ) are 
the renormalized normal and anomalous Green's func- 
tions, respectively. The resonance peak appears only in 
limited regions of the q-space around (jr, it) and (jr, —ir). 
We show in Fig. 8 the temperature variation of the q- 
integrated intensity spectrum. It is desirable to have the 
corresponding neutron inelastic scattering experiments 
on single crystals in future although such an experiment, 
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Fig. 7. Dispersion and broadening of the resonance peak or the 
spin excitons. (a) T/t = 0.0070, (b) T/t = 0.0143. 



with a deuterized sample, does not seem to be very easy. 

The nuclear spin-lattice relaxation rate,_can be cal- 
culated from the following general formulaEa' using the 
FLEX dynamical susceptibilities: 

,2r 



1 



Q 



N, 



A, 



o 



,2 Imx + (<?,^o) 

U) 



(4) 



where 7^ is the nuclear gyro-magnetic ratio, A q the 
Fourier ^-component of the hypcrfine coupling constant, 
and lo the resonance frequency. We may neglect the 
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Fig. 8. g-integrated spectrum of the imaginary part of the dy- 
namical susceptibility. 
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Fig. 9. Calculated nuclear spin-lattice relaxation rate 1/Ti (log- 
arithmic plots). The dashed line indicates a T 3 behavior. Inset 
shows 1/T\T (normal plots). 



q-dependence of the hyperfine coupling constant which 
is anisotropic at the site of 13 C. According to ref. 8 
the principal values parallel and perpendicular to the di- 
rection of the 13 C= 13 C bond are a + 2B and a — B, 
respectively, with a = 1.3kOc//iB and B = 2.1kOe//iB 
for one spin per molecule; thus these values should be 
divided by 2 for one spin per dimer. For experimental re- 
sults under the magnetic field parallel to the conducting 
layer we may approximately replace \A q | in eq. (4) with 
[(a + 2B) 2 + (a- B) 2 ]/8, since the C=C bond direction is 
nearly perpendicular to the layer. The result of calcula- 
tion using this value for the coupling constant is shown in 
Fig. 9. The T 3 dependence of the relaxation rateJpelpw 
T c as observed experimentally is well reproduced£mE£P 
The order of magnitude seems also reasonable, in view of 
the approximate nature of the calculation. The observed 
values of 1/TiT above T c has a peak at T = T* ~ 50K. 
Above this temperature T* the system shows insulating 
behavior while the system is metallic below T*. Thus 
above T* we had better use the Heisenberg model with 
the Anderson superexchange. 



For T* > T > T c the observed values for \jT{T de- 
crease with decreasing temperature, reminiscent of the 
pseudo-spin gap behaviors in high-T c cuprates. This be- 
havior is hard to understand within the present model 
and approximation. One possible explanation may be 
given by considering a coexistence of spin density wave 
(SDW) fluctuations and a charge density wave (CDW) 
or its slow fluctuations, although further studies are nec- 
essary before drawing any conclusion. 

In summary, we have studied the spin fluctuation- 
induced superconducting state of quasi-two dimensional 
k-(BEDT-TTF) 2 A@ compounds by using a half-filled 
triangular Hubbard model within the FLEX approxima- 
tion, extending the previous study above T c . The energy 
gap of (x 2 — y 2 )-type develops below T c more rapidly 
than in the BCS model. Resonance peaks in the dynam- 
ical susceptibility are predicted in limited regions of the 
q-space around (tt, tt) and (ir, —%). The observed T 3 be- 
havior of the nuclear spin-lattice relaxation rate 1 /T± is 
well reproduced with a reasonable order of magnitude. 
We have also predicted the doping concentration depen- 
dence of T c in these compounds. Calculations of the SC 
and AF instability points in the U/t vs. t'/t plane show 
that the former is favorable only for t'/t larger than 
~ 0.3, indicating that an antiferromagnetic instability 
dose not always accompany the d-wave superconductiv- 
ity. 

Finally we would like to emphasize that the present 
mechanism seems to be the only available mechanism 
to describe the anisotropic superconductivity in_£[uasi- 
two dimensional K-type BEDT-TTF compounds.^ The- 
oretical results seem to be generally consistent with 
available experimental results except for the pseudo-spin 
gap behavior of 1/Ti above T c )l3 Further experimental 
and theoretical investigations are desired to confirm this 
mechanism. 

We would like to thank Dr. S. Nakamura and Dr. T. 
Takimoto for helpful discussions. 
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